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Abstract 

We test the method of analytic continuation from imaginary to real chemical potential in two- 
color QCD, which is free from the sign problem. In particular, we consider the analytic contin- 
uation of the critical line to real values of the chemical potential. 
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1. Introduction 

Understanding the phase diagram of QCD on the temperature-chemical potential (T, /j,) 
has many important implications in cosmology, in astrophysics and in the phenomenology 
of heavy ion collisions. It is known that the fermion determinant in QCD becomes complex 
in presence of a non-zero chemical potential, thus preventing us from performing standard 
Monte Carlo simulations. Various techniques have been exploited to circumvent this 
problem. Among these techniques we consider here the use of an imaginary chemical 
potential for analytic continuation [1-16]. It is very important to have some control on 
the method of analytic continuation and a direct test of this method in simpler models 
can be very useful. To this purpose we explored the case of 2-color QCD where, for any 
value of the chemical potential, the fermion matrix is real. In ref. [9] we found, by looking 
at the behavior in both real and imaginary chemical potential of the Polyakov loop, the 
chiral condensate and the fermion density, that the method gives reliable results, within 
appropriate ranges of the chemical potential, and that a considerable improvement can 
be achieved if ratio of polynomials are used to interpolate data with imaginary chemical 
potential [9,10] . 
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Fig. 1. Critical couplings obtained from the chiral susceptibility, together with a linear fit in (a^t) 2 to 
data with fj? < 0. The datum at /3 C = 1.30 (diamond) is taken from Ref. [14]. 

On the other hand one of the most important applications of analytic continuation is 
the determination of the critical line or of critical surfaces for small values of /i [1-6,11,12]. 
The theoretical basis in this case is not as straightforward as for physical observables 
and relies on the assumption that susceptibilities, whose peaks signal the presence of 
the transition, be analytic functions of the parameters on a finite volume [1,2]. Direct 
tests of the method are even more important in this case, therefore we extended [17] our 
analysis of two-color QCD to the study of the (pseudo)critical line. As in usual QCD 
simulations, we will determine locations of the critical line for /i 2 < and interpolate 
them by suitable functions to be continued to /i 2 > 0. The prediction obtained at real /i 
will then be compared with direct determinations of the transition line. 

2. Lattice simulations 

Numerical simulations have been performed using hybrid Monte Carlo with dt = 0.01 . 
Typical statistics collected corresponds to 20k trajectories. Simulations have been done 
using resources of the INFN apeNEXT computing center. 

Our aims are finding the best interterpolation of critical line at imaginary chemical 
potential and testing the analytic continuation of the critical line from imaginary to real 
chemical potential. We considered SU(2) gauge theory with Nf — 8 staggered fermions 
with mass am = 0.07 on a 16 3 x 4 lattice [17]. In Fig. 1 we present results for the 
extrapolation of the critical line to /i 2 > compared with the direct determination of 
the critical coupling (3 C . As one can see there is a discrepancy between the extrapolated 
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Fig. 2. Our preliminary results for the critical line in SU(3) at finite isospin on a 8 3 X 4 lattice and 
am = 0.1, together with a linear fit in (a/i) 2 to data with ft 2 < 0. 

critical line and the direct determination of /3 c (/i 2 ). One can infer that or the critical 
line is not analytic in the whole range of /x 2 considered here or that the interpolation at 
/i 2 < is not accurate enough to correctly reproduce the behavior at /i 2 > 0. A common 
fit to all available data with a third order polynomial in /j.2 suggests the second possibility 
may be the correct one. In order to go deep inside this question we have started new 
simulations using a different mass (am = 0.2) on a 16 3 x 4 lattice and mass am = 0.07 
on a 16 3 x 6 lattice. We have also considered the test of the analytical continuation of 
the critical line the case of SU(3) with finite isospin and Nf = 8 staggered fcrmions with 
mass am = 0.1 on a 8 3 x 4 lattice. Fig. 2 displays our preliminary results for the analytic 
continuation of the critical line in the case of SU(3) at finite isospin density. As one can 
see the discrepancy between the extrapolated critical line and the direct determination 
of PdlJ 2 ) is less severe. 



3. Conclusions 



We have tested in two-color QCD the analytic continuation of the critical line in 
the (T, /i) plane from imaginary to real chemical potential. We have found that the 
critical line around /j, = can be described by an analytic function. Indeed, a third order 
polynomial in /i 2 fits all the available data for the critical coupling. However, when trying 
to infer the behavior of the critical line at real n from the extrapolation of its behavior 
at imaginary /z, a very large precision would be needed to get the correct result. In the 
case of polynomial interpolations there is a clear indication that high-order terms play 
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a relevant role at /i 2 > but are less visible at /x 2 < 0, this calling for an accurate 
knowledge of the critical line in all the first Roberge- Weiss sector. This scenario could 
be peculiar of two-color QCD. If confirmed in other theories free of the sign problem, 
such as QCD at finite isospin density, then one should seriously reconsider the analytic 
continuation of the critical line in the physically relevant case of QCD at finite baryon 
density. 

For the reason above we have started a systematic program of investigations both in 
two-color QCD and in SU(3) QCD at finite isospin. Our preliminary results suggest that 
it may be useful to revisit the interpolations used in QCD (first order polynomial in /jr) 
and make the effort of determining accurately at least one more term in the polynomial 
fit. 
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